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Different classes of physical systems with sizeable electron-phonon coupling and lattice distor- 
tions present anomalous resistivity behaviors versus temperature. We study a molecular lattice 
Hamiltonian in which polaronic charge carriers interact with non linear potentials provided by local 
atomic fluctuations between two equilibrium sites. A path integral model is developed to select 
the class of atomic oscillations which mainly contributes to the partition function and the electrical 
resistivity is computed in a number of representative cases. We argue that the common origin of 
the observed resistivity anomalies lies in the time retarded nature of the polaronic interactions in 
the local structural instabilities. 
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I.INTRODUCTION 

While a sizeable electron-phonon coupling is prerequi- 
site for polaron formation, it is not generally true that a 
system with strong electron-phonon interactions (g > 1) 
has to display polaronic features. In a simple metal with 
wide bare bands and long range forces like lead, for in- 
stance, intersite off-diagonal hopping terms prevent trap- 
ping of the charge carriers and onset of a polaronic state. 
In strong coupling systems with a narrow electronic band 
D the antiadiabatic inequality D/TiQ < 1 (cD being the 
characteristic phonon frequency) is likely fulfilled thus 
ensuring the stability of the unit comprising electron and 
phonon cloud that is, the small polaron |lj . Also in adi- 
abatic regimes small polarons can form although a real 
space broadening of the quasiparticle has to be expected, 
at least in three dimensional systems. As a consequence, 
the adiabatic polaron mass renormalization is less strong 
than that occuring for the antiadiabatic polaron once the 
strength of the electron-phonon coupling has been fixed 
H . Theoretical investigation on these issues has become 
intense after that a polaronic mechanism has been pro- 
posed ID as a viable possibility to explain high T c super- 
conductivity with its peculiar transport properties. Evi- 
dence has been recently provided Q for the polaronic na- 
ture of carriers also in perovskite manganites with colos- 
sal magnetoresistance (CMR) and a strong coupling of 
carriers to Jahn- Teller lattice distortions has been sug- 
gested H to explain the resistivity peak located at the 
Curie temperature |(| . Although polarons are invoked to 
account for different resistivity behaviors in systems with 
some lattice distortions |7j, the common origin, if any 
does exist, of such behaviors remains unexplained and a 
unifying theory is lacking. In general, a strong electron- 
phonon coupling implies violation of the Migdal theorem 
H and polaron collapse of the electron band with the ap- 
pearance of time retarded interactions in the system. The 
path integral method p|JlO|] has proved successful in deal- 
ing with this problem since a retarded potential naturally 
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emerges in the exact integral action. Both the Frohlich 
p] , ^2| and the Holstein [|l3|-[l5| polaron properties have 
been studied by using path integrals techniques. 

In this paper we propose a path integral approach to 
the problem of a polaron scattered by a local lattice insta- 
bility arising from an anharmonic phonon mode. Small 
polarons are assumed to exist by virtue of the strong 
overall electron-phonon coupling and independently of 
the single anharmonic mode. As an example, the lat- 
ter could be associated with the motion of apical oxy- 
gen atoms in Y Ba-iCu^O-i superconductor where a c- 
axis polarized high frequency phonon couples to in-plane 
charge carriers On the other hand, the highly cor- 
related motion of in-plane holes and c-axis lattice vibra- 
tions could in itself give origin to a small polaron [fj"7| . 
We model the local instability by a double well poten- 
tial in its two state configuration, a Two Level System 
(TLS), and the one dimensional atomic path between the 
two equilibrium minima is taken as time dependent [fl8| . 
Retardation effects are thus naturally introduced in the 
full partition function and the path integral method per- 
mits to derive the effective time (temperature) dependent 
coupling strengths which control the resistivity. 



II. THE PATH INTEGRAL MODEL 

Our starting point is contained in the following Hamil- 
tonian 

H {r) = e( 5 )St( T )2( r ) + ^ w q 4(r)a q (r) + H tls {t) 

q 

( AQ(r) 
- \XQ{t) 

flint (r) = -2\Q(t)(*(t)c(t) 
Q(t)= -Q + —(r-t t ). (0.1) 

To 



(h tls {i 



t is the time which scales as an inverse temperature 
according to the Matsubara Green's function formalism. 
Hq (t) is the free Hamiltonian made of: a) a polaron cre- 
ated (distroyed) by (t)(c(t)) in an energy band e(g) 
whose width decreases exponentially by increasing the 
strength of the overall electron-phonon coupling constant 
g, e(g) = Dexp(—g 2 ); b) a diatomic molecular lattice 
whose phonon frequencies w q are obtained analytically 
through a force constant approach; c) a Two Level Sys- 
tem in its symmetric ground state configuration due to an 
anharmonic phonon mode. Q(r) is the one dimensional 
space-time hopping path followed by the atom which 
moves between two equlibrium positions located at ±Q D . 
t is the bare hopping time between the two minima of 
the TLS and U is the instant at which the ith-hop takes 
place. One atomic path is characterized by the number 
2n of hops, by the set of t, (0 < i < 2n) and by r Q . We 
have implicitly assumed that the class of r-linear paths 
yields the main contribution to the full partition function 
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of the interacting system [[19| . The closure condition on 
the path is given by: (2n — l)r s + 2utq = f3, where j3 is 
the inverse temperature and r s is the time one atom is 
sitting in a well. The interaction is described by Hint (t) 
with A being the coupling strength between TLS and po- 
laron, \Q(r) is the renormalized (versus time) tunneling 
energy which allows one to introduce the r dependence 
in the interacting Hamiltonian pof . 

Following a method previously developed |nj in the 
study of the Kondo problem, we multiply XQ(t) by a 
fictitious coupling constant s (0 < s < 1) and, by dif- 
ferentiating with respect to s, one derives the one path 
contribution to the partition function of the system 



Z{n,U 
Z 



= -2A f ds [ dTQ(r) lim G(t,t') s 



(0.2) 

where, Zq is the partition function related to Hq and 
G(r, t ) s is the full propagator for polarons satisfying a 
Dyson's equation: 

G(t,t') s = G°(r,r) + s [ dyG (T,y)\Q(y)G(y,r' ) 

(0.3) 



The polaronic free propagator G° can be derived ex- 
actly in the model displayed by eqs.(l). We get the full 
partition function of the system by integrating over the 
times U and summing over all possible even number of 
hops: 



71=0 



f 13 dt 2n 


f t2 - T0 dti 




■ / — exp 


Jo T o 


Jo T a 



Z T = Z o y^ y / — — • • / — -exp —pE(n,U,To) 



0E(n, t l ,r a )=L-(K A + K R )J2{ ^ j («-4) 

with E(n,ti,To) being the one path atomic energy. L, 
which is a function of the input parameters, is not essen- 
tial here while the second addendum in eq.(4) is not local 
in time as a result of the retarded polaronic interactions 
between successive atomic hops in the double well po- 
tential. K A and K R are the one path coupling strengths 
containing the physics of the interacting system. K A 
(negative) describes the polaron-polaron attraction me- 
diated by the local instability and K R (positive) is re- 
lated to to the repulsive scattering of the polaron by 
the TLS. Computation of E(n,ti,To) and its derivative 
with respect to To shows that the largest contribution to 
the partition function is given by the atomic path with 
t s = 0. The atom moving back and forth in the dou- 
ble well minimizes its energy if it takes the path with 
the highest To value allowed by the boundary condition, 
that is with (To)maa: = (^uKbT)^ 1 . This result, which 
is general, provides a criterion to determine the set of 
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dominant paths for the atom at any temperature. Then, 
the effective interaction strengths < K A > and < K R > 
can be obtained as a function of T by summing over n 
the dominant paths contributions: 



<K A >= -(\Q ) 2 B 2 exp(2j2A< l )Y, A ^lf 

q q 

< K R >= -f3(XQ ) 3 B 3 exp(3 £ A n ) £ A^j\f (0.5) 

q q 

with: B = {np{l) — l)exp(—g^2 tl ctgh(f3ui cl /2)) and 
A q = 2 5 ^7V q (7V q + i). 

Nq is the phonon occupation factor and np(e(g)) is the 

Fermi distribution for polarons. / = J2n=i( T o)max an d 
N is the cutoff on the number of hops in a path. The 
particular form of {To) max suggests that many hops paths 
are the relevant excitations at low temperatures whereas 
paths with a low number of hops provide the largest con- 
tribution to the partition function at high temperatures. 
Since the effective couplings which determine the resis- 
tivity depend on (ro)^ aa ., hence on iV~ 3 (through /), a 
relatively small cutoff (N ~ 4) ensures numerical conver- 
gence of eqs.(5) in the whole temperature range. On the 
other hand the non retarded term L in eq.(4) has a slower 
1/N behavior, therefore a larger cutoff should be taken 
at low temperatures where computation of equilibrium 
properties such as specific heat is strongly influenced by 
many hops atomic paths between the minima of the dou- 
ble well excitations. 

Our lattice Hamiltonian is made of diatomic sites 
whose intramolecular vibrations can favor trapping of 
the charge carriers [f3| . The miramolecular frequency 
wo largely influences the size of the lattice distortion as- 



sociated with polaron formation 21 while the dispersive 
features of the phonon spectrum controlled by the in- 
termolecular couplings are essential to compute the po- 
laron properties both in the ground state and at finite 
temperatures. The range of the intermolecular forces is 
extended to the second neighbors shell since these cou- 
plings remove the phonon modes degeneracy (with re- 
spect to dimensionality) at the corners of the Brillouin 
zone thus permitting to estimate with accuracy the rel- 
evant contributions of high symmetry points to the mo- 
mentum space summations. 

In a simple cubic lattice model with first and second 
neighbors molecular sites interacting via a force constants 
pair potential, the w q are: 



2 a + 37 + 4(5 1 rz 

<= M + MV" 2+F ^ 

^x,y,z — r ) / ^(.CxCy ~f~ Sx$y H~ Qc^z H~ ^x^z ~t~ CyC z 

+s y s z )) + 2a-/(c x + c y + c z ) + AaS(c x c y + c y c z ) 
+2j5(3c y + 2c z + c x + 2cos(q x - q y - q z ) 

+2 S Z ) + C x C 2 y + S x S 2y + CyC 2z + S y S 2z ) 

+2<5 2 (2 + 2c x c z + 2c z (c x c 2y + s x s 2y ) + c 2y + c 2z ) (0.6) 
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where, c x = cosq x , s x = sinq x , C2 y — cos2q y etc. 
The intramolecular force constant a is related to loq 
by Wq = 2a /M with M being the reduced molecular 
mass. 7 and S are the intermolecular first and sec- 
ond neighbors force constants, respectively. Let's define 
oj 2 = 7/M and u>2 = S/M. Then the characteristic fre- 
quency Co, which we choose as the zone center frequency, 
is: Co 2 = ui 2 + zuj 2 + z nnn uj' 2 , z is the coordination number 
and z nnn is the next nearest neighbors number. 

We turn now to compute the electrical resistivity due 
to the polaronic charge carrier scattering by the impurity 
potential with internal degree of freedom provided by the 
TLS. Assuming s wave scattering, one gets [p2|,p| 

p = posin 2 r] 

po = — 5- (0.7) 

ne 2 v 2 F (N /V) Ti 

where, n s is the density of TLS which act as scatterers, 
vf is the Fermi velocity, V is the cell volume, iVo is the 
electron density of states, e is the charge and Ti is the 
Planck constant. The phase shift ij to the electronic wave 
function at the Fermi surface is related to the effective 
interaction strengths < K A > and < K R >. 

The input parameters of the model are six that is, the 
three molecular force constants, g, D and the bare en- 
ergy A<5o- Qo can be chosen as ~ 0.05 A consistently 
with reported values in the literature on TLS's which are 
known to exist in glassy systems p3| , p4[ , amorphous met- 
als 25 1, A15 compounds || and likely in some cuprate 



superconductors |L7 26|. In these systems the origin of 
the TLS's is not magnetic. The bare electronic band 
D is fixed at O.leV. In Fig.l, we take a rather large Co 
setting the system in moderately adiabatic conditions 
and an intermediate {g — 1) to very strong {g > 3) 
electron-phonon coupling regime. The broad resistivity 
maximum developing at low temperatures in the interme- 
diate regime clearly signals that the TLS's are at work 
here |27|] while, by increasing T, the hopping time short- 
ens and the incoming polarons cannot distinguish any 
more the TLS internal degree of freedom, then diago- 
nal scattering prevails and a quasi linear resistivity be- 
havior emerges at T > 200i^. At larger g the polaron 
becomes progressively smaller in real space and the oc- 
curence of the self trapping event (at g ~ 3) is marked 
by an abrupt increase in the effective mass, in any di- 
mensionality |2^,^]. As a consequence, the resistivity 
becomes non-metallic and its absolute values are strongly 
enhanced. Let's fix g = 1 which ensures polaron mobil- 
ity and tune (see Fig. 2) the TLS-polaron coupling A. We 
see that a resistivity peak located at T ~ 150-ftT arises 
at A > 700mel/A _1 with height and width of the peak 
being strongly dependent on A hence, on the TLS energy. 
The low T resistivity still displays the maximum at the 
unitary limit while the high T (T > 300!^) behavior can 
be metallic like ( A < SOOmeVA" 1 ) or semiconducting 
like (A > mOmeVA- 1 ). A ~ 700 - SOOmeVA^ 1 cor- 
responds to a TLS energy of ~ 35 — AQmeV which is 
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comparable to the value of the bare polaron energy band 
e(g = 1) ~ 37meV. In this picture, the resistivity peak 
has a structural origin and it can be ascribed to resonant 
TLS-polaron scattering with effective attractive and re- 
pulsive interaction strengths becoming of the same order 
of magnitude at T ~ 15QK. While the peak does not shift 
substantially by varying the strength of the intermolec- 
ular forces the height of the peak turns out to be rather 
sensitive to those strength hence, to the size of the quasi- 
particle. As an example, by doubling w\ with respect to 
the value in Fig. 2 (at fixed loq and w 2 ), the normalized 
resistivity peak (A = lOOOmeWl -1 curve) drops to 0.65 
while by taking u)\ — 60meV the peak disappears. The 
effect of the second neighbors force constant is relatively 
less strong: by doubling U2 (at fixed luo and u)\) with 
respect to fig. 2, the resistivity maximum becomes 0.76 
while u>2 = 30meV yields a peak value of 0.66. 

At larger A's the resonance peak is higher and broader 
since an increasing number of incoming polarons can be 
off diagonally scattered by the TLS. However, the ap- 
pearance of this many body effect mediated by the local 
potential does not change the position of the peak either, 
which instead can be shifted towards lower temperatures 
by reducing substantially ujq hence the Q value. In Fig. 3, 
we set luq = 50meV. As previously discussed Q the 
intermolecular coupling lo\ should be ~ wo/2 to ensure 
a correct application of the molecular lattice model with 
first neighbors interacting sites. Here we are extending 
the range of the forces to the second neighbors shell ac- 
cordingly, lj\ is slightly reduced to the value 20meV while 
0J2 is switched on and fixed at lOmeV. This choice of pa- 
rameters allows us: i) to treat correctly the ground state 
polaron properties versus dimensionality pq |, ii) to set 
the characteristic phonon frequency of the present three 
dimensional lattice model at u> ~ 75meV which is the 
energy of c-axis polarized phonons due to apical oxygen 
vibrations coupled to the holes in the Cu-0 planes of 
Y Ba-zCuzO-i-?, |]l6| . Although our simple cubic lattice 
does not account for the details of the structural effects 
of YBCO we are in the appropriate range of parame- 
ters to capture the main features of the lattice polarons 
scattered by local instabilities in those compounds. 

Anharmonic features of the oxygen modes have been 
recognized to be larger in underdoped samples |3C| ] and 
doping dependent polaron formation |n| has been cor- 
related to distortions of the oxygen environment. Our 
model predicts (Fig. 3) a large resisitivity peak whose 
height is strongly increased with respect to Fig. 2 (com- 
pare the A = lOOmeVA^ 1 cases) due to the fact that 
the polaron effective mass is heavier when lower energy 
phonons build up the quasiparticle. The non metallic 
behavior at T larger than ~ 90K reminds of the anoma- 
lous c-axis resistivity observed in underdoped high T c 
superconductors p2|j . We believe that the semiconduct- 
ing like p c in underdoped high T c superconductors can 
be ascribed to anharmonic potentials due to oxygen dis- 
placements strongly coupled via A (as in Fig. 3) to po- 
laronic carriers. To attempt a comparison with experi- 
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ments we need to fix the residual resistivity through the 
ground state parameters given in eq.(7). By extrapolat- 
ing to T = the normal state data on Y BaiCu-sOi-?, 
one derives p(T — 0) ~ 3mf2cm, hence the experimental 
peak value p ~ 20mQcm observed in the 80K compound 
Y Ba<iCuiP§.%7 [B2| can be reproduced in my model by 
A ~ 1300meV/L which corresponds to a local mode en- 
ergy of ~ 65meV, in fair agreement with the measured 
energies of phonons strongly coupled to the charge carri- 
ers. 

In the range 80 A < T < 200 A the A ~ UOOmeVA- 1 
curve fits rather well the data of ref. §§]. At T = 200K, 
the experimental p is ~ lOmflcm and the calculated 
value is ~ 8mf2cra. In the range 200A < T < 300A, the 
experimental p is essentially flat while my results still ex- 
hibit a dp/ dT < behavior. In general, the resistivities 
here obtained have an exponential temperature depen- 
dence (at T larger than the resonance) consistent with a 
polaronic hopping motion whereas the experimental data 
show a variety of sample dependent non metallic behav- 
iors for the out of plane resistivity |}3| in underdoped 
compounds. 

A resistivity peak at T ~ 250A" is observed in CMR 
materials || and a splitting of the Mn states due to a 
local lattice distortion has been proposed as a possible 
mechanism. In our description, only extended polarons 
dragging a cloud of high frequency optical phonons (see 
Fig. 4) can produce a peak at T ~ 200 — 250A which 
is however broader than the experimental one. High 
phonon frequencies imply reduced mass renormalization 
and, consistently, weak overall electron-phonon coupling 
values. Extended polarons move therefore in wide en- 
ergy bands and the resonance effect can take place only 
if the TLS splitting energy is large enough. Infact, by 
taking D = Q.heV and A ~ bOOOmeVA- 1 (Fig.4), we get 
the peak at g < 1 whereas, in a small polaron regime 
(g > 1) the resonance is lost and the peak is suppressed. 
Although more specific models accounting for magnetic 
field effects are required for the CMR systems, we argue 
that the strong coupling of polarons to structural distor- 
tions is likely relevant to their transport properties. 

III. CONCLUSION 

The path integral method provides a powerful tool to 
describe non local ( in time ) scattering of the charge car- 
riers by local (in space) potentials due to structural dis- 
tortions. Our charge carriers are lattice polarons which 
exist in the diatomic molecular lattice by virtue of a size- 
able overall electron-phonon coupling. The structural 
instability is due to an anharmonic atomic motion be- 
tween two equilibrium positions which typically could be 
~ Q.1A apart. I have derived the full partition function 
of the interacting system and obtained the temperature 
dependent effective couplings by summing over the class 
of energetically favored atomic hopping paths. We sug- 
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gest that the retarded nature of the interactions is the key 
ingredient to explain some anomalous transport proper- 
ties observed in real materials. A resistivity peak is ob- 
tained as a consequence of resonant scattering of polarons 
strongly coupled to local double well potentials. 



FIG. 1. Electrical Resistivity normalized to the resid- 
ual (T = 0) resistivity for five values of electron-phonon 
coupling g. The bare TLS energy AQo is 6.5meV. The 
force constants which control the phonon spectrum are: 
u)o — lOOme V, uii = 20meV, u>2 = WmeV 



FIG. 2. Electrical Resistivity normalized to the residual re- 
sistivity for five values of the polaron-TLS coupling A. g = 1. 
The intra- and intermolecular force constants are as in Fig.l 
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FIG. 3. Electrical Resistivity normalized to the residual 
resisitivity for nine values of the polaron-TLS coupling A. 
g = 1. o>o = 50meV, u>i = 20meV,u>2 = WmeV 

. The experimental data (X) are taken from ref. |Q . 



FIG. 4. Electrical Resistivity normalized to the residual 
resistivity for four overall couplings g. A = 5000meVA~ . 
ujq = WOmeV, lui = 80meV, uj 2 = 20meV 
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